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A theoretical model of electromigrative, diffusive and convective transport polymer- 
gcl composites is presented. Bulk properties are derived from the standard clcctroki- 
netic model with an impenetrable charged sphere embedded in an electrolyte-saturated 
Brinkman medium. Because the microstructurc can be carefully controlled, these mate- 
rials arc promising candidates for enhanced gcl-clectrophoresis, chemical sensing, drug 
delivery, and microfluidic pumping technologies. The methodology provides solutions for 
situations where perturbations from equilibrium arc induced by gradients of electrostatic 
potential, concentration and pressure. While the volume fraction of the inclusions should 
be small, Maxwell's well-known theory of conduction suggests that the model may also 
be accurate at moderate volume fractions. In this work, the theory is used to compute 
ion fluxes, electrical current density, and convective flow driven by an electric field ap- 
plied to an homogeneous composite. The electric-field-induced (electroosmotic) flow is 
a sensitive indicator of the inclusion ^-potential and size, electrolyte concentration, and 
Darcy permeability of the gel, while the electrical conductivity is most often indepen- 
dent of the polymer gel and is relatively insensitive to characteristics of the inclusions 
and electrolyte. 



1. Introduction 

Gel-electrophoresis is widely used to sort macromolecules based on their size and elec- 
trical charge. Selectivity to size is achieved by adjusting the permeability of the gel 
(e.g., agarose or polyacrylamide) through the concentrations of monomer, cross-linker, 
catalyst and initiator used in the gel synthesis. Molecular sorting based on other char- 
acteristics, such as receptor-ligand binding affinity, requires the gel to exhibit specific 
physicochcmical activity. One way to achieve this in a controlled manner is to embed 
surface- functionalizcd particles (e.g., biological cells, and synthetic polymer or silica 
spheres) in a conventional polymer gel. Accordingly, this work seeks to quantify the 
influence of surface charge on transport driven by gradients of chemical and electrostatic 
potential. While the task is simplified to some extent by limiting the analysis to sim- 
ple electrolytes whose mobilities are unhindered by the polymer gel, the methodology 
provides a significant step toward a theory that also accounts for hindered transport of 
larger electrolyte ions (e.g., proteins and DNA fragments). 

This work also provides a quantitative interpretation of novel diagnostic tests — analogous 
to well-established microelectrophoresis and conductivity measurements — that probe the 
surface charge or ^-potential of immobilized colloids in electrolytes where the particles 



2 



would otherwise aggregate. Attractive particle-interaction potentials arise when the solu- 
tion pH approaches the isoelectric point of the parti cle-electrolyte inter face, or the surface 
charge is sufficiently well screened by added salt ijRussel et al\\l98^i . By immobilizing 
colloids in an (ideally) inert (uncharged) polymer gel, at a pH and ionic strength where 
the interactions arc repulsive, the pH and ionic strength may be varied without inducing 
coagulation. 

Membranes of sintered glass beads (without intervening polymer) have long been used 
in ion -selective electrodes, and, more recently, as electroosmotic pumps (e.g.. lYao et all 
Their simple design (no moving parts) and high-pressure low-flow characteristics 
are ideally suited to microfluidic applications. Filling the void space with a permeable, 
uncharged polymer gel, as proposed in this work, will increase viscous dissipation and, 
therefore, diminish pumping efficiency. Nevertheless, because applications are envisioned 
where poor pumping efficiency might be tolerated in view of other attributes, a quanti- 
tative analysis of electroosmotic pumping is undertaken here. 

The ch arge on particles di spersed in an electrolyte endows them with an electrophoretic 
mobi lity (|Russel et a/ll989|) . Th eoretical inte rpretation of the mobility JO'Brien fc White] 
19781). se dimentation potential l)Savillelll982(l . low-frequen cy conductivity (ISavillelll979t 
Q'Brienlll98ll) . dielectric response (complex c onductivity) foeLacev fc Whitdll9'8l|) , and 



electroacoustic response (dynamic mobility) l)0'Brienlll988l Il99d|) is widely used to in- 
fer the surface charge and, therefore, to indicate dispersion stability. Closely related are 
streaming-potential and streaming-current devices, which are used to infer the charge 
on macroscopi c substrates, and the charge dens ity and permeability of porous plugs and 
coatings (e.Q.. lHunte^l200ltlDukhin et oi.ll2004h . 

This paper addresses a new but related problem in which impenetrable spheres with 
surface charge are randomly dispersed and immobilized in a permeable, electrolyte- 
saturated polymer gel (see figure^]). The transport processes that take place with the ap- 
plication of average (macroscale) gradients of electrostatic potential, electrolyte concen- 
tration and pressure are derived. While transport of electrolyte ions is relatively straight- 
forward to calculate in an homogeneous, uncharged gel, charged inclusions disturb the 
applied fields, and applied fields disturb the equilibrium state of the diffuse double lay- 
ers, so the resulting fluxes reflect a complex coupling of electromigration, diffusion, and 
convective transport. 

Electrokinetic th eories are often based on the standard electrokinetic model l(Overbeekl 
119431 lBoothlll95fi) . whereby continuum equations governing the electric field, mobile 
charge (microions), mass, and momentum are solved with appropriate boundary con- 
ditions. The principal difficulty usually lies in capturing double-layer polarization and 
relaxation at surfaces whose radius of curvature is comparable to or smaller than the 
equilibrium double-layer thickness (Debyc length). For 'bare' and 'soft' (polymer-coated) 
parti cles, polarization and relaxation can be addressed with novel numerical methodolo- 
gies (lO'Brien fc Whitelll978tfHill et al.ll2003l) . In this work, a numerically exact solution 
of the problem is achieved for the low-volume-fraction limit where particle interactions 
can be neglected. The methodology also neglects quadratic and higher-order perturba- 
tions to the equilibrium state, providing asymptotic coefficients that characterize the 
far-field (power-law) decays of the velocity disturbance and perturbations to the equilib- 
rium electric field and ion concentrations. In turn, the asymptotic coefficients arc linked 
to bulk properties of the composite. The analysis resembles Maxwell's theory for the 
effective conductivity of a dilute, random configuration of spherical inclusions. 

Dipolar disturbances^ arise in the limit where the inclusions arc uncharged. Then, with 



f Here, dipolar refers to axisymmetric disturbances, without a source, that satisfy Laplace's 
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the application of a uniform electric field, the perturbed electrostatic potential ip' reflects 
the non-conducting (impenetrable) surface of the inclusions, satisfying V 2- 0' = 0. Simi- 
larly, when subjected to a bulk concentration gradient, the perturbed ion concentrations 
n' satisfy V 2 n' = 0. In both cases, there is no convective flux, because the electrolyte is 
everywhere electrically neutral. With a uniform pressure gradient (Vp), the problem sim- 
plifies to the flow of an incompressible Newtonian fluid through a Brinkman medium with 
impenetrable inclusions. The (solenoidal) velocity u satisfies ?yV 2 M — Vp — (rj/£ 2 )u = 
0, where p is the pressure, r\ is the fluid viscosity, and £ is the Brinkman screening 
length (square root of the Darcy permeability). When the inclusion radius a ^> i, 
u = — (£ 2 /?/)Vp with V 2 p = (Darcy flow), and the drag force on the inclusions 
is 2irqU a? £~ 2 , where U = — (£ 2 /if)(^7p) + 0(4>) is the average fluid velocity. In this 
limit, the inclusion contribution to the average drag force (per unit volume) of the com- 
posite is n2irr]Ua 3 £~ 2 = (3/2)(j)riU£~ 2 1 where n is the inclusion number density and 
<j> = n(4/3)vra 3 is the volume fraction. With the same far-field velocity, the drag force on 
each inclusion is clearly much greater than the Stokes drag force QirriaU . Note also that 
the velocity disturbances u' = u — U = — (l/2)a Uj(6ijT — 3xiXjr~ 5 ) decay as r~ 3 
when the distance from the inclusions r > a > /. 

When the inclusions are charged, a diffuse layer of mobile counter-ions envelops each 
inclusion, and electroneutrality demands that the net charge in the diffuse layers balances 
the immobile surface charge. With an applied electric field, the electrical body force 
within the diffuse layers drives an 'inner' flow which, in turn, drives an 0(<f>) 'outer' flow, 
with an O(l) contribution due to an imposed pressure gradient. The electric-ficld-induced 
velocity disturbances 'pump' fluid through the polymer gel without exerting a net force 
on the composite. While the velocities are very low, composites with an a 1 cm 2 cross- 
section can produce velocities of several microns per second in a microchannel. Note that 
a pressure gradient is necessary to overcome the drag required to pump fluid through 
an external network. However, the pressure-driven contribution to the flow is often small 
compared to the electric-field-induced flow. Under these conditions, the electric-field- 
induced flow rate is practically independent of the pressure gradient, and the maximum 
pressure gradient that can be sustained is limited by the strength of the composite and 
its support. 

Ion fluxes manifest in an electrical current and, hence, a measurable electrical conduc- 
tivity. Follo wing earlier treatments of th e low-frequency conductivity of dilute colloidal 
dispersions llSa,vil1elh979l lO'Brienlll Q8lL the incremental contribution of the inclusions 
to the effective conductivity is calculated. These results link conductivity measurements 
to the particle surface charge density, for example. Because the fluxes are dominated by 
elcctromigration, the conductivity is not significantly influenced by the gel. 

The model also provides the effective diffusivities of electrolyte ions when the bulk 
electrolyte concentration varies slowly in space and time. In a forthcoming paper (to be 
published elsewhere), two important situations are examined: (i) bulk diffusion in the 
absence of an average electric field, and (ii) bulk diffusion with an electric field yielding 
zero current density. The former provides a simple setting in which to demonstrate the 
influence of the inclusions on the effective ion diffusion coefficients, whereas the latter pro- 
vides the particle contribution to the membrane diffusion potential, which is well known 
in the fields of membrane biology, electrochemistry, and electrochemical engineering. In 
both cases, the particle contribution to these bulk properties may be comparable to or 
larger than in the absence of inclusions. 

equation; these take the form otjdr^ 1 / dxj — —ajXjr~ 3 , where r = \x\ and Oi are the compo- 
nents of a constant vector. 
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The paper is organized as follows. We begin in with a description of the electroki- 
netic model. First, the composite microstructure and an apparatus for comparing theory 
and experiment are described. The subsections therein present the electrokinetic trans- 
port equations ( fc!2.1[l . which are used to compute linearly independent solutions of the 
single-particle (microscale) problem f H2.2[) . Asymptotic coefficients from solutions of the 
microscale problem capture far-field decays of perturbations to the equilibrium state. 
These are used to calculate bulk ion fluxes in ^2.31 and to derive an average momentum 
equation in M2.4I Results are presented in S|3]for composites with negatively charged inclu- 
sions in polymer gels saturated with NaCl electrolyte. The subsections therein examine 
the incremental pore mobility f H3.1|) . electroosmotic pumping fi j3.2|) . incremental pres- 
sure gradient and, finally, species fluxes ( H3.4JI and electrical conductivity f HH.5[) . 
A brief summary follows in 21 



2. Theory 

The microstructure of the composites considered in this work is depicted in figure ^ 
The continuous phase is a porous medium comprised of an electrically neutral, electrolyte- 
saturated polymer gclf. Polyacrylamide gels are routinely used for the electrophorctic 
separation of DNA segments in aqueous media. Their porosity may be controlled by 
adjusting the average densities and ratio of the monomer (acrylamide) and cross-linker. 
In this work, the hydrodynamic permeability is characterized by the Darcy permeability 
£ 2 (square of the Brinkman screening length), which reflects the hydrodynamic size a s 
and concentration n s of the polymer segments. In turn, these reflect the degree of cross- 
linking and the affinity of the polymer for the solvent. 

Embedded in the polymer are randomly dispersed spherical inclusions. In model sys- 
tems, the inclusions are envisioned to be monodispcrsc silica or polymeric spheres, which 
typically have radii in the range a = 10 nm-10 /im and bear a surface charge when 
dispersed in aqueous media. The surface charge density may vary with the bulk ionic 
strength and pH of the electrolyte. In this work, however, the surface charge is to be 
inferred from the bulk ionic strength and surface potential Because the inclusions are 
impenetrable with zero surface capacitance and conductivity, no-flux and no-slip bound- 
ary conditions apply at their surfaces. 

Note that the mobile ions whose charge is opposite to the surface-bound immobile 
charge are referred to as counter-ions, with the other species referred to as co-ions. For 
simplicity, the counter-charges, i.e., the dissociated counter-ions, are assumed indistin- 
guishable from the electrolyte counter-ions. Surrounding each inclusion is a diffuse layer 
of mobile charge, with Debye thickness k _1 and excess of counter-ions. As described be- 
low, the layer structure is calculated from the well-known Poisson-Boltzmann equation. 

In this work, the polymer gel is assumed not to hinder ion motion. For larger ions and 
dense polymer gels, the influence of the network on the diffusive, electromigrative and 
convective fluxes may be modeled with an equation of motion for an ion 

y(u-v) -j'v + f = 0, (2.1) 

where the first term represents the hydrodynamic drag due to relative motion; u and 
v are the (average) fluid and ion velocities, and 7 is the friction coefficient. The second 
term approximates the force exerted by the polymer gel on the ion, with the friction 
coefficient 7' reflecting the relative size of the ion and polymer interstices. The third 

f A gel refers to a network of polymer chains that is cross-linked so as to exhibit a solid-like 
(elastic) response to an applied stress. 
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Figure 1. Schematic of the microscale system under consideration. Charged, impenetrable 
inclusions [solid circles) with radius a ~ 10 nm-10 /im are embedded in a continuous polymer gel 
(solid filaments) saturated with aqueous electrolyte. Diffuse double layers (dash-dotted circles) 
with thickness kT 1 ~ 1-100 nm are perturbed by the application of an average electric field 
— (VV>), pressure gradient (Vp), or electrolyte concentration gradient (V%). The Brinkman 
screening length £ ~ 1-10 nm that specifies the Darcy permeability £ 2 of the gel is often small 
compared to the radius of the inclusions. 

term accounts for electrical and thermal (Brownian) forces, depending on the time-scale 
of interest. The unhindered ion velocity is v° = u + f /j, so the hindered velocity may 
be written v = v°j/(j + 7'). Therefore, under steady conditions, the ion conservation 
equation V- (nv) = is independent of 7' (n is the ion number density). When u = 0, for 
example, the ion will diffuse (or clcctromigrate) with an effective diffusivity (or mobility) 
D e = 1)7/(7 + 7') 1 where D is the unhindered diffusivity. It follows that v = v°D e /D, 
so the hindered flux equals the unhindered flux multiplied by the ratio of the hindered 
to unhindered ion diffusivitics (mobilities). 

Now consider the influence of hindered ion migration on the fluid momentum conser- 
vation equation. From Eqn. (|2.1|l . the (hydrodynamic drag) force exerted by an ion on 
the solvent is 

7(«-«) = (/-y«)7/(7 + 7'), (2-2) 
where / = — zeV^ is the electrical force on the ion (z is the valance, e is the elementary 
charge, and ip is the electrostatic potential). It follows that the net force (per unit volume) 
exerted by the ions on the fluid is 

N 

- njWjU + zjeV^/Dj, (2.3) 

where the sum is over all N ion species. Clearly, as 7^/7; — * when the hindrance of 
the polymer is negligible, ions transfer their electrical force to the fluid. As Jj/^fj — ► 00, 
however, immobilized ions transfer the electrical force to the polymer, so the net force 
exerted by each ion on the fluid becomes —jjU. 

A simple apparatus to realize the conditions under which the theory may be applied 
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Figure 2. Schematic of the macroscale system under consideration. A polymer gel embedded 
with spherical charged inclusions (see figure^ separates (by length L) two reservoirs containing 
electrolyte with different species concentrations (rij°(z = 0) and n°°(z — L)) and, possibly, 
pressures (p(z = 0) and p(z = L)). Electrodes on each side of the 'bridge' permit an electric 
field to be applied and the differential electrostatic potential Atf) = ^(z = L) — tp(z = 0) to be 
measured. The walls of the bridge are impenetrable and non-conducting. 



is depicted in figure [3 The composite bridges two reservoirs, each, in general, with a 
different electrolyte concentration and pressure. Electrodes are placed at each end of the 
bridge, so either a uniform electric field can be established or an average electric field 
strength measured. The channel is to realize constant ion fluxes under steady or quasi- 
steady conditions. This makes the averaged microscalc transport equations easier to solve, 
but, in general, the averaged equations apply to macroscale fluxes in two- and three- 
dimensional geometries when the (average) inclusion number density, ion concentrations, 
electric field and fluid velocity vary slowly in space and time. 

2.1. The electrokinetic transport equations 

The transport equations and boundary conditions arc presented here in dimensional form. 
They comprise the non-linear Poisson-Boltzmann equation 

N 

e e s V 2 V = ~J2( n i ~ n j) z i e > ( 2 - 4 ) 
i=i 

where e and e s are the permittivity of a vacuum and dielectric constant of the electrolyte; 
nj are the concentrations of the jth mobile ions with valences Zj\ and tp and e are the 
electrostatic potential and elementary charge. In this work, the polymer is uncharged, so 
the fixed charge density is zerof . 

Transport of the mobile ions is governed by 

6miaj{u - Vj) - ZjeVip - kTVlnnj = (j = 1, N), (2.5) 

where aj are Stokes radii of the ions, obtained from limiting conductances or diffusivitics; 
rj is the electrolyte viscosity; u and Vj are the fluid and ion velocities; and kT is the 
thermal energy. 

Ion diffusion coefficients, which arc adopted throughout this paper, are 

Dj = kT/(6nr)aj). (2.6) 

t For convenience, the valence of the fixed charge is set opposite to that of its respective 
(mobile) counter-ion in Eqn. I|2.4fl : this defines the concentration nj. 
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As usual, the double-layer thickness (Debye length) 

k- 1 = ^/kTe s e /(2Ie 2 ) (2.7) 
emerges from Eqns. (|2.4|) and (|2.5|) where, 

N 

7=(l/2)]>»r (2.8) 

i=i 

is the bulk (average) ionic strength, with n°° the bulk ion concentrations. 
Ion conservation demands 

drij/dt = - V • ( nj Vj) (j = l,...,N), (2.9) 

where t is the time, and the ion fluxes 

are obtained from Eqn. Ij2.5jl . 

Similarly, momentum and mass conservation require 

N 

p s du/dt = nV 2 u - Vp- {r)/£ 2 )u - ]P n^-eV^ (2.11) 

3=1 

and 

V-w = 0, (2.12) 

where p s and u are the electrolyte density and velocity, and p is the pressure. Note 
that — {rj/£ 2 )u represents the hydrodynamic drag force exerted by the polymer on the 
electrolyte. The Darcy permeability I 2 (square of the Brinkman screening length) of the 
gel may be expressed as 

f = l/[n s (r)67ra s F s ] = 2a 2 J[9 ( f> s (r)F s (0 s )}, (2.13) 

where n s (r) is the concentration of Stokes resistance centers, with a s and F s ((j) s ) the 
Stokes radius and drag coefficient of the polymer segments. In this work, n s is constant, 
but, in general, may vary with radial position r from the center of each inclusion. Note 
also that the Brinkman screening length is adjusted according to Eqn. I|2.13[) by varying 
the (uniform) polymer segment density with St okes radius «.<, = 1 A. T he drag coefficient 
F s is obtained from a correlation developed bv lKoch fc Sanganl |l999j) for random fixed 
beds of spheres. While the microstructurc of a polymer gel is clearly not the same as that 
of a random bed of spheres, the model is intended to capture the significant influence of 
hydrodynamic interactions on the permeability when the volume fraction of polymer is 
not small. In this work, however, only the reported values of £ are relevant. For example, 
the value £ ~ 0.951 nm, which is adopted for the principal set of results tabulated 
below, reflects a polymer segment concentration n s that yields I = 1 nm according to 
Eqn. (|2.13|l when the Stokes radius a s ~ 1 A and the drag coefficient F s = 1. Because 
the hydrodynamic volume fraction 4> s = n s (4/3)7ra^ > 0, F s ((f> s ) > 1 and, hence, £ is 
slightly less than the targeted value. 

2.1.1. Inner (particle surface) boundary conditions 

Either the equilibrium surface potential £ or surface charge density a may be spec- 
ified. Because the surface (r = a) is assumed impenetrable with zero capacitance and 
conductivity, the surface charge is constant, permitting no- flux boundary conditions for 



each (mobile) ion species. As usual, the no-slip boundary condition applies. It follows 
that (inner) boundary conditions are either 

V> = C at r = a (2.14) 

or 

e s e Vip\ 0Ut ■ h - e p e Vip\ in ■ n = -er at r = a, (2-15) 

with 

rijVj ■ n = at r = a (2.16) 

and 

u = at r = a. (2.17) 
where h = e r is an outward unit normal and e p is the particle dielectric constant. 

2.1.2. Outer (far- field) boundary conditions 

Neglecting particle interactions requires far-field boundary conditions 

ip ^ -E ■ r &s r ^ oo, (2.18) 

rij — * + Bj ■ r as r — ► oo, (2-19) 

and 

u -> U as r -> oo, (2.20) 

where £7, £?j and f/ are, respectively, a constant electric field, constant species concen- 
tration gradients, and constant far-field velocity. 

2.2. Solution of the equations 

2.2.1. Equilibrium state 

When E = Bj = U = 0, equilibrium is specified according to 

N 

e e s V 2 V° - - - nfce, (2.21) 
j'=i 

= V • [DjVn$ + z.e^nyr] (2.22) 

and 

JV 

= - Vp° - ^ n^j-eV^ , (2.23) 

with boundary conditions 

^° = C at r = a, (2.24) 

e s e o V-0°|out • e r - e p e o Vip \ in ■ e r = -a at r = o, (2.25) 

n°[L>., Vn° + z,-e-£n?VV°] ■ e r = at r = a, (2.26) 

and 

V>° -> as r -► oo, (2.27) 
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n° -> n°° as r -► oo. (2.28) 

2.2.2. Linearized perturbed state 

Perturbations to the equilibrium state (above) are introduced via 

$ = <ijp - E ■ r + %l>\ (2.29) 

rij • B ■ r ■ „',. (2.30) 

and 

p = p" + P-r+p\ (2.31) 

where the first terms on the right-hand sides denote the equilibrium values, with the 
primed quantities denoting perturbations. Note that P is the far-field pressure gradient 
required to sustain a far-field velocity U = — (£ 2 /n)P. 
With 'forcing' 

X = Xe z , (2.32) 

where X £ {E, Bj, C/}f , the linearized perturbations are symmetric about the z-axis 
(9 = 0) of a spherical polar coordinate system, taking the forms 

i' = t/j(r)X ■ e r (2.33) 



n'j = nj(r)X ■ e r (2.34) 

and 

u = U + u', (2.35) 

wherej 

u' = V x V x h(r)X 

= -2{h r /r)(X ■ e r )e r - (h rr + h r /r)(X ■ e e )e e . (2.36) 

The perturbations satisfy 

N 



3=1 



and 
where 



V-j,=0, (2.38) 



= -^-(VnJ + Bj) - Zje^in'j + B 3 ■ r)VV>° 

-Zje^W -E) + n°j(U + u'), (2.39) 

t In general, a linear combination of these variables. 
£ Equation 12.361 guarantees a solenoidal (incompressible) velocity field, which permits the 
momentum equation to be solved by applying the curl V x , thereby eliminating the pressure 
and yielding a scalar equation for the single non-zero component of the vorticity V x u = w</,e</,. 
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and 

TV 

77V V - Vp' - {v/e 2 ){U + u') - J^nJjZj-efV^' - E) 

3=1 

N 

~ UK" + B o ■ ^ZjeVtp = (2.40) 
V-w' = 0, (2.41) 

with boundary conditions 

e s e (VV/ - E)\out ■ e r - e p e (Vip' - E)\ in ■ e r = at r = a, (2.42) 



D, 



and 



+z j e T ^n i (Vi(j' -E)-n°JU + «')] ■e r = 0atr = a ! (2.43) 

u' = U at r = o, (2.44) 

V/ -> (X ■ e r )D x /r 2 as r -> 00, (2.45) 

n^-> (X ■ e r )Cf /r 2 as r -> 00, (2.46) 

u' -2(C x /r 3 )(X • e r )e r - (C x /r 3 )(X • e e )e e as r -» 00. (2.47) 

In the far field, the velocity disturbance u' is proportional to the gradient of p' , which, 
like the electrostatic potential and ion concentrations, is dipolar. Accordingly, u' decays 
as r~ 3 , and D x and C x will often be referred to as the strength of the electrostatic and 
concentration polarization (or dipole moments) induced by the field X £ {E, Bj, U}. 

The dimensions of the asymptotic coefficients D x , C x and C , which depend on 
the respective X £ {E, Bj, U}, are easily worked out by inspecting Eqns. I|2.45(l - (|2.47|l . 
For convenience, dimcnsionless values are presented in the tables below with a, u* = 
e s e (kT/e) 2 /(rja), 21 and kT/e as the scales for length, velocity, ion concentrations, and 
electrostatic potential, respectively. 

2.2.3. Superposition 

The equations are solved using a numerical methodology developed hv lHill et al\ J2003|) 
for the electrophoretic mobility of polymer-coated colloids. Solutions with E, Bj and U 
set to arbitrary values can be computed, provided Ylj=i z j^j = to ensure an electrically 
neutral far-field. However, when vV species are assembled into M electroneutral groups 
(e.g., electrolytes or neutral tracers), each with far- field gradient Bk (k = 1,...,M), it 
is expedient to compute solutions with only one non-zero yalue of E, Bk or P. T hen, 
arbitrary solutions can be obtained by linear superposition l|0'Brien fc Whitelll978l) . 

An index k' is required to identify the (electroneutral) group to which the jth species 
under consideration is assigned. Careful consideration of the electrolyte composition and 
ion valences is required to ensure consistency. For z-z electrolytes it is convenient to 
set Bj = Bk, whereas for a single 2-1 electrolyte (e.g., CaC^) is it satisfactory to set 
Bj = B^ /\zj\. For the relatively simple situations considered in this work, the (single) 
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electroneutral group is NaCl, so M = 1 with k = k' = 1, and j = 1 and 2 for Na + and 
Cl~, respectively. 

Note that C Bk , for example, is the asymptotic coefficient for the perturbed concen- 
tration of the jth species induced by the kth. concentration gradient B^, whereas G Bk 
(without a subscript) denotes the asymptotic coefficient for the flow induced by B^. 

For neutral species, the concentration disturbance produced by a single impenetra- 
ble sphere yields Cf k ' = (l/2)a 3 , otherwise Cf k = (k ^ k'). Clearly, the asymptotic 
coefficients for charged species, whose concentration perturbations are influenced by elec- 
tromigration, are not the same as for neutral species; for ions, C^ k ' — ► (l/2)a 3 as \C\ — > 0, 
however. 

With co-linear forcing and bulk electroneutrality, linear superposition gives far-field 
decays 

M 

ip' -> (l/r 2 )[ED E + V B k D Bk + UD u ](e z ■ e r ) as r -> oo, (2.48) 
fc=i 

M 

{l/r 2 )[ECf + Y. BkC f" + UC fK e * ■ e r) as r ^ oo, (2.49) 

k=l 

M 

v! -> -(2/r 3 )[^C B + + l/C^Ke, • e r )e r 

k=l 

M 

-(l/r 3 )[EC E + B kC Bk + UC u ]{e z ■ e g )e g as r -» oo. (2.50) 

fe=i 

This work is pr imarily concerned with situations where only one X S {E, Bj, U} is 
applied. Following rO'Brien fc White! l|l978|) . these are referred to as the (E), (B) and (U) 
(microscale) problems. Algebraic or differential relationships between the averaged fields 
may be applied to ensure zero average current density, for example. The next section 
relates (microscale) E, Bj and U to the averaged (macroscale) fields, e.g., — (Vi/>}, 
(Vn 3 -) and (it), in dilute composites. 

2.3. Averaged (bulk) fluxes 
Here we calculate the average flux of the jth species 

(jj) = V- 1 J j s dV, (2.51) 

where the volume of integration includes the continuous and discrete phases. If the size of 
the representative elementary volume is between the micro- and macro-scales, the result 
is equivalent to sampling the flux (at a point) over all micro-structural configurations 

(ensemble a verage) . 

Following ISavillel (Il970h and [o] Brienl l)l98l[) . the averaging can be accomplished by 
adding and subtracting the flux 

-DjVrij - Zje^nfVip + n°°u (2.52) 
kT J J 

from the integrand in Eqn. H2.51J) . This yields the macroscopic electromigrativc, diffusive 
and convective fluxes in the absence of inclusions, plus an integral whose integrand is 



and 
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exponentially small beyond the diffuse double layers, i.e., 

(J'i) = -Zje^nfiyt) - DjiVnj) + nf{u) 



-V 



1 J ^ e ^r n T v ^ + D * Vn ^ - n ? u + hW- ( 2 - 53 ) 



Applying the divergence theoremf and noting that V • jj = 0, the volume integral in 
Eqn. (|2.53|) becomes 

/ z i e ^ n T^ dA ~ J Zje^rijiV^ ■ n}rdA 
+ j DjUjhdA - J DjlVrij ■ n]rdA 

+ J (rij - nf)(u ■ n)rdA, (2.54) 

where the surface integrals enclose the inclusions and their respective equilibrium double 
layer, with n directed outward, into the fluid. 

For dilute composites, i.e., when n(4/3)7r(a + k^ 1 ) 3 -C 1, the integral over a represen- 
tative volume V equals nV integrals with a single particle at r = 0, each with h = e r . 
Therefore, noting that n® — n°° is exponentially small as r — > oo, and that n'j and Bj ■ r 
are odd functions of position, the fluxes become 



(j 3 .) « nf (u) - z ^nf{Vi>) Dj(V nj ) 
+nZje Ri n °° j _ (w' ■ r)]e r dA 

+nDj / [n'j - (Vn'j ■ r)]e r dA. (2.55) 

J r — >oc 

Superposing solutions of the independent single-particle problems, the microscale elec- 
tromigrative and diffusive contributions to the average flux, i.e., the last two terms in 
Eqn. (|2~55|) . arej, 



D J oo 

nzje—rij 



f [ip' - (W • r)]e r dA 

J r— >oo 



' r — >oo 

M 



n^z 3 e^-nf[ED E + V B k D B * +UD U ], (2.56) 



kT 

fc=i 



and 



. M 

nDj / [nj - (Vn'j ■ r)]e r dA = nAirDj[ECf + ^ B k Cf k + UCf}. (2.57) 

Jr^oo fe=1 



t Note also that J ad!/ = f. cca • ndA — f v iV • adV, where a represents an arbitrary 
vector field. 

X Note that r(a ■ e r ) has been written as (a ■ r)e r because r = re r , where a represents an 
arbitrary vector field. 
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Substituting these into Eqn. I|2.55|l gives 

« nf(u) - z 3 e^nf(V^) D^Vn,) 

d j 

— -7 

kT 

+n^D j [ECf + J2 B kCf k + UCf\. (2.58) 



n. M 
+n^z e-^nf[ED E + £ B k D B « + UD U ] 



k=l 
M 



k=l 



Note that the average fluxes are now expressed in terms of the asymptotic coefficients 
from at most 2 + M independent single-particle problems, each of which is solved 'exactly' 
in this work. 

2.4. Averaged (bulk) momentum conservation equations 

In general, an average velocity (u) is produced by the application of an average pres- 
sure gradient (Vp), electric field — (V^>), or concentration gradients (Vrij). This section 
relates these to the asymptotic coefficients emerging from the single-particle problem. 
Averaging the fluid momentum equation gives (see Appendix 0) 

« -<Vp) - (r)/£ 2 )(u) + V W 2 (u) - (pV^j) - n(f d ) (2.59) 

where p = p° + p' is the charge density and (f d ) is the average (hydrodynamic) force 
exerted by the fluid on the inclusions. Note that inertia is neglected, as arc hydrodynamic 
and electrostatic interactions; the analysis is therefore limited to small volume fractions 
n(4/3) 7 r(a + K- 1 ) 3 < 1. 

Similarly to the average fluxes, let us adopt the single-particle problem to evaluate 
(f d ). For a single inclusion in an unbounded polymer gel, 



f [-(P-r+p')5 + 2r)e\ ■ e r dA 

J r—a 

= [ [-(P • r +p')S + 2r]e] ■ e r dA - f [(r)/£ 2 )u + pVi/j]dV, (2.60) 

J r — >oo J r — a 

where e = (1/2)[Vm+(Vu) t ] and 5 is the identify tensor. Since Vw = 0, u(r = a) = 0, 
and Vp = —(ri/£ 2 )u and u 1 ~ r~ 3 as r — > oo, Eqn. (|2.6U|) can be written 

(f d ) « - f \p'- (Vp' ■ r)]e r dA - f pV^dV. (2.61) 

J r — >oo J r—a 

Beyond the double layer, 

p' = -(v/t 2 ) (2/r' 3 )C x (X-e r )dr> = -(l/r 2 )( V /f)C x (X ■ e r ), (2.62) 



and because the integrand of the volume integral in Eqn. 12.61f) is exponentially small 
there, 



(pVip) w n I pVipdV. 

J r — a 



(2.63) 



Therefore, 

M 

(f d ) « { v /£ 2 )4ir[EC E + B k C Bk + UC U ] - n^ipVip) (2.64) 



fc=i 
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and, hence, Eqn. Ij2.59(l becomes 



M 



P3 -(Vp) - {ri/£ 2 )(u) + nV 2 {u) - n(r]/£ 2 )4n[EC E + B k C Bk + UC U }. (2.65) 

k=l 

Note that, in addition to (f d ), an electrical force (/ e ) and a mechanical-contact force 
(f m ) ac t on the inclusions. Accordingly, static equilibrium requires 

en = -<r> - (f d ) 

M 

« -{ti/1 2 )A-k[EC e + B kC Bk + UC U }. (2.66) 

k=l 

In the absence of charge, for example, — (f m ) is equal to the drag force on a sphere 
embedded in a B rinkman medium with viscosity r\ and Darcy permeability £ 2 . Indeed, 
iBrinkmanl (|l947^ solved this problem exactly, obtaining 

(f rn ) = -2^0(0/ £) 2 [l + 3(£/a) + 3{£/a) 2 ] (( = 0, (j> -> 0), (2.67) 

which shows that 

2C u /a 3 = 1 + 3{£/a) + 3{£/a) 2 (C = 0,0^0). (2.68) 

Note that the drag force approaches the Stokes drag QnrjaU as £/a — > oo. When £/a — > 0, 
however, the drag approaches 27ra 3 (j) / £ 2 )U because the surface traction is dominated by 
the pressure (dipole) arising from the outer Darcy flow: V 2 p' = with u = -(l 2 /rj)Vp'. 

2.5. Averaged (bulk) equations for unidirectional transport 
With all average fluxes in the z-direction, mass and momentum conservation require 
constant (it) and, hence, 

M 

(Vp) = -(v/i 2 )(u) - c/>(3/a 3 )( V /£ 2 )[EC E + ]T£ fc C B " + UC U }. (2.69) 

fe=i 

Similarly, the (steady) average species conservation equations V • (jj) = require con- 
stant average fluxes 

(jj) = nf(u) - Zj e^nf(Vi>) - Dj(Vnj) 

n. M 
+<j>(3/a 3 )Zje^n?[ED E + £ B fc D B * + UD U ] 



k=l 
M 



kT 

+cf>(3/a 3 )D J {EC E + J2 B k Cf* + UCf]. (2.70) 

fc=i 

Note that V • (V?A) = in an electrically neutral composite with uniform dielectric 
permittivity, so the average electric field is also constant. 

The averages can be expanded as power series in the inclusion volume fraction e.g., 
(u) — > Uq + (f>Ui + 0(4> 2 ). Therefore, since the microscale equations (asymptotic coeffi- 
cients) are accurate to 0((j)), the notation is condensed by writing, for example, (u) = U , 
where it is understood that U = U + (j)Ui + O((j) 2 ). Clearly, E, Bj and U in Eqns. 12.6911 
and i|2.70|) need only to include the 0(1) contribution to their respective average, e.g., 
U — > Uq. The following notation is adopted for the other averaged quantities: Jj = (jj), 
P = (Vp), Bj = (Vnj), E = -<VV). 
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With one electrolyte (M = 1) and, recall, bulk electroneutrality, there are TV + 4 
independent variables (E, U, P, Bk (k = 1), Jj (j = 1, TV)) with TV + 1 independent 
equations (see Eqns. I|2.(j9[) and ^2J[)\i ). Clearly, three independent variables must be 
specified for a unique solution. 

For clarity, the results presented below involve a 1-1 electrolyte (NaCl), mostly with 
only one non-zero forcing variable. It is important to note that, because the equations 
are linear, solutions for any combination of forcing variables may be constructed. For 
example, a forthcoming paper establishes the electric field strength required to maintain 
a constant electrolyte flux — driven by a bulk concentration gradient across a membrane — 
with zero bulk current density. 

3. Response to an electric field 

Results are now presented for the application of an electric field in the absence of 
average pressure and concentration gradients. These conditions prevail when measuring 
the electrical conductivity, for example, and they provide a relatively simple setting 
in which to study the influence of inclusions on bulk electrokinetic transport. Steady 
homogeneous conditions are assumed, neglecting the influence of electrode polarization 
and electrochemical reactions. Accordingly, the average velocity from Eqn. I|2.69[l is 

U = -<t>{i/a 3 )EC E + O(0 2 ), (3.1) 

and the average ion fluxes from Eqn. H2.70|) are 

Jj = z ^ n T E + 4>^/o a )z J e^nfED E 

+<f>{3/a 3 )D J EC E + (t>{i/a 3 )nfEC E + 0((j) 2 ). (3.2) 

Asymptotic coefficients arc provided in tabled for a composite with Brinkman screen- 
ing length £ sa 0.951 nm and inclusion radius a = 100 nm; the ^-potentials and (three) 
ionic strengths span experimentally accessible ranges. With a positive electric field (E > 
0), the counter-ions (Na + ) migrate toward the 'front' of the inclusions, inducing a pos- 
itive electrostatic dipole moment D E > 0. When the ^-potential is low, however, the 
dipole moment reflects the dielectric polarization required to maintain an impenetra- 
ble interface, so the dipole strength approaches the Maxwell value D E = — (l/2)a 3 (for 
non-conducting spheres) as |£| — > 0. The positive concentration dipole moments C E > 
reflect the combined influences of electromigration, diffusion, and electroneutrality. As 
expected from Eqn. Ij2.47|l . C E < 0, because the electrical force on the fluid and, hence, 
the resulting electroosmotic flow are forward (U > 0). 

3.1. Incremental pore velocity 
As suggested by Eqn. I|3.1[) . the ratio 

U/{E<t>) = -3C E /a 3 , (3.3) 

which is termed the incremental pore mobility, provides a convenient measure of the 
electroosmotic pumping capacity. When multiplied by the electric field strength and 
particle volume fraction, the values in the last column of tabled f° r example, yield the 
0(cj)) average velocity that prevails in the absence of an applied pressure gradient. This 
section examines how the strength of the flow is related to the ^-potential and size of the 
inclusions, the ionic strength of the electrolyte, and the permeability of the gel. We will 
see that the pore mobility is significantly influenced by polarization and relaxation, so 



Table 1. Dimensionless asymptotic coefficients (see Eqn. (13. 2H 1 and incremental pore mobility (Eqn. 13,31 1 for bulk electromigration of NaCl 
in a Brinkman medium with charged spherical inclusions: a = 100 nm; £ ~ 0.951 nm; T = 25°C; Di » 1.33 x 10" 9 m 2 s- 1 (Na+); 
D 2 « 2.03 x 10" 9 m 2 s _1 (CP); u* = e s e {kT / ef / (no) « 5.15 x 10" 3 m s _1 . 



Ce/(fcT) 


D E /a 3 


CfkT/(2Ia 3 e) 


C E kT/(u*a 4 e) 


U/{E<j>) = -3C E /a 3 






(7 = 1,2) 




((nm S _1 )/(V cm- 1 )) 


/■£(2 = 1 


j n oc w in-6 tv^I 1 — 1 
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— o.oZ X iu 
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+1.96 x 10° 
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I = 9.25 x 10~ 4 mol l" 1 
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-4.74 x 10" 1 


7.30 x 10~ 2 


-1.20 x lO" 3 
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I = 9.25 x 10~ 2 mol r 1 








-1 
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the qualitative form of the relationship (with a given gel permeability) is similar to the 
classical electrophoretic mobility of dispersions ([O'Brien fc Whitall973) . 

Consider the case in tablenwith na = 100 and ( = —lkT/e. With E = 2 V cm -1 and 
= 10~ 2 , the pore velocity is U ~ 0.84 nm s _1 , which is clearly very slow. If, however, 
this flow is directed from a composite with a 1 cm 2 cross-section into a microfluidic 
channel, then it is not unreasonable to amplify the velocity by four orders of magnitude, 
yielding a (modest) average velocity of 8.4 /mi s _1 . In this example, the permeability (£ 2 = 
0.951 2 nm 2 ) is relatively low, so higher velocities may be achieved with a comparable 
(weak) electric field and inclusion volume fraction. Note that a stronger electric field 
between (platinum) electrodes separated by a distance L ~ 5 mm, say, produces hydrogen 
and oxygen bubbles. In conventional electroosmotic pumps, much higher el ectric field 



2003). For the purpose of accurately determining the pore mobility, however, higher 
electric field strengths are, perhaps, unnecessary. 

The pore mobility is shown in the left panels of figure[3]as a function of the ^-potential 
for various values of na, with inclusion radii a = 10, 100 and 1000 nm (top-to-bottom 
panels), an d Brinkman screening len gth £ ks 3.11 nm. Similarly to the electrophoretic mo- 
bility (e.g.. lO'Brien fc Whitdll978h . at low to moderate ^-potentials the pore mobility 
provides a one-to-one connection between the (measured) pore velocity and the surface 
charge. Mobility maxima arise from polarization (by electromigration) and relaxation 
(by diffusion) of the equilibrium double layer. As suggested by earlier theoretical stud- 
ies examining the role of polarization and relaxation on the electrop horetic mobility of 
polymer-coated particles (|Savillell200f)HHiljl2004tlHiil fc Savillell2005r) . these calculations 
clearly demonstrate that polarization is driven by electromigration, since convection is 
extremely weak when the particles are immobilized in a polymer gel. Note that theoreti- 
cal studies of electroosmotic flow in micro-porous membranes do not reveal such maxima, 
because the u nderlying microscale m odel comprises (effectively) straight channels with 
charged walls ftao fc Santiagoll2003h . 

As expected, the (incremental) pore mobility tends to increase with ^-potential at fixed 
ionic strength, and increases with ionic strength at fixed ^-potential. Both trends reflect 
the increasing charge required to maintain a constant surface potential when varying the 
ionic strength. For colloids whose surface charge is independent of electrolyte concentra- 
tion, the ^-potential increases with decreasing na. Accurate semi-empirical expressions for 
this relations hip (obtained from solutions of the Poisson-Boltzmann equation) are read- 
ily available IjRussel et aZ.lll989j) . In general, however, the dependence of surface charge 
density on ionic strength and pH is exceedingly difficult to predict, and must therefore 
be determined empirically for specific interfaces (e.g., see lYao et al 1l2003i for silica in the 
presence of KC1). 

To draw a closer connection to experiments, the pore mobility is shown in the right 
panels of figure |31 as a function of the ionic strength with three constant surface charge 
densities spanning two orders of magnitude. Note that the inclusion radii are the same 
as in the corresponding left panels. Because the surface charge is fixed, the ^-potential 
(dashed lines, right axis) decreases with increasing ionic strength, but the particle size 
does not significantly influence the £-potential. Because the average velocity reflects the 
combined influence of all particles in the composite, U is expected to be proportional 
to the O(noo? ~ acfi/a) (average) counter-charge density. Therefore, balancing the cor- 
responding 0(Eacj)/a) electrical force with the 0(rjU/£ 2 ) Darcy drag force gives a pore 
mobility U/(<f>E) ~ a£ 2 /(r]a). Indeed, comparing the mobility axes (left sides) of the 
panels on the right-hand side of figure |21 indicates that the mobility is, at least approxi- 
mately, inversely proportional to the inclusion radius. Again, at low ionic strengths, when 
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the ^-potential is high, polarization and relaxation significantly complicate this simple 
interpretation. 

To highlight the influence of polymer-gel permeability, the pore mobility is shown in 
figure 0] as a function of the Brinkman screening length £ for various ^-potentials. With 
a particle radius a = 100 nm, the ionic strength / sa 0.0925 mol 1 _1 yields na = 100 and 
nl > 1 at most values of £. Now the mobility increases as £ m with exponent m{£) in the 
range 1-2, indicating that viscous stresses and Darcy drag balance the electrical body 
force. Note that the mobility increases linearly with the ^-potential when |£| is small, 
and, again, mobility maxima are evident when |£| « 6kT/e. Pore mobilities are shown in 
figureElat a much lower ionic strength / sa 9.25 x 10~ 6 mol l^ 1 yielding na = 1. Now, as 
expected, the mobility increases linearly with £ 2 , since k£ < 1. The monotonic increase 
with ^-potential is because the surface charge densities are low and, hence, polarization 
is weak. 

3.2. Electroosmotic pumping 

Recall, the (incremental) pore mobility was defined with zero average pressure gradi- 
ent, and therefore neglects the pressure differential Ap = PL required to pump fluid 
through an external network. This section briefly addresses coupling the composite and 
electrodes — which together are referred to as an electroosmotic pump — to a microflu- 
idic network. The analysis briefly considers the force exerted on the composite, electrical 
power consumption, and pump efficiency. 

Let us consider a closed loop, where fluid in the composite exits from one side, passes 
through a microfluidic network, and returns to the other side. To ensure that all electrical 
current flows through the composite, we assume that the electrical resistance of the 
network is much greater than that of the pump, which is realized when the length (area) 
of the external network L e (d 2 ) is much greater (smaller) than that of the composite L 
(A). Next, assuming laminar viscous flow, the pressure-drop through the network may 
be written 

Ap » V c(L e /d 4 )Q, (3.4) 

where d 2 is the (characteristic) channel cross-sectional area, Q is the volumetric flow 
rate, and c is an 0(1) constant that reflects the shape and length the network sections. 
Equating the network pressure drop to the pump characteristic emerging from Eqn. 1)2. 65J) . 

Ap/L = -(v/^U - 0{3/a 3 )( V /£ 2 ){EC E + UC U ), (3.5) 

gives 

[rjc(L e /d 4 )(A/L) + (r)/£ 2 ) + <t>(3/a 3 )(r,/£ 2 )C u ]Q/A = -^/a^/l^EC* , (3.6) 

where A and L are the composite cross-sectional area and length. 

The contribution of the back-flow, as represented by the asymptotic coefficient C u 
in Eqn. (|3.6|l . is obtained from the (U) problem. This and other asymptotic coefficients 
are provided in table El for a composite with Brinkman screening length £ « 0.951 nm 
and inclusion radius a = 100 nm; again, the ^-potentials and (three) ionic strengths 
span experimentally accessible ranges. Here, the average force exerted by the polymer 
on the inclusions (f m ) = ~4:Tt(t]/£ 2 )C u U (Eqn. I|2.66|l ) is independent of the surface 
charge. Note that C u = (l/2)a 3 when the only contribution to the force is due to Darcy 
flow (£Ja <C 1); the (constant) value C u /a 3 w 0.514 in table |3 which reflects a small 
viscous contribution (£/a = 0.00951), is precisely the value given by Brinkman's theory 
(Eqn. H2.67|l ). The flow-induced electrical and concentration polarization, as represented 
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Figure 3. The incremental pore mobility U / (Ecf>) with inclusion radii a = 10 (top), 100 (middle) 
and 1000 nm (bottom): aqueous NaCl at T = 25°C; I ~ 3.11 nm. Left panels show the mobility 
as a function of the (scaled) ("-potential C,e/(kT) for various (scaled) reciprocal double-layer 
thicknesses ko, = 1, 10, 30, 100, 300 and 1000. Right panels show the mobility (solid lines, left 
axis) and (-potential (dashed lines, right axis) as a function of the ionic strength with constant 
surface charge densities —a ~ 0.179, 1.79 and 17.9 fj,C cm -2 (increasing upward at high ionic 
strength) . 
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Figure 4. The incremental pore mobility U/(E<f>) as a function of the Brinkman screening 
length £ for various (scaled) (-potentials — (e/(fcT) = 1, 2, 3, 6 (solid lines) 7, 10 (dashed 
lines): aqueous NaCl at T = 25°C; a = 100 nm; na = 100 (I « 0.0925 mol l" 1 ). An electric 
field is applied in the absence of average pressure and concentration gradients. The maximum 
velocity is achieved when \Qe/(kT) ~ 6. 




Brinkman permeability (nm ) 

Figure 5. The incremental pore mobility U/(E<j>) as a function of the Darcy permeability I 2 
for various (scaled) ("-potentials — (e/(fcT) = 1, 2, 3, 6 (solid lines) 10 (dashed line): aqueous 
NaCl at T = 25°C; a = 100 nm; ko = 1 (I « 9.25 x 10 -6 mol l" 1 ). 

by D y and Cj 7 , are related to the streaming potential and streaming current, and are 
included here only for future reference. 

When the dominant resistance comes from the composite itself (the second term on 
the left-hand-side of Eqn. I|3.6|l ). the pump performance curve simplifies to 

Q pa -(j>(3/a 3 )AEC E , (3.7) 

which is clearly independent of the applied load. Neglecting constraints imposed by elec- 
trolysis, for example, the maximum length of the composite may be set by consideration 
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Table 2. Dimensionless asymptotic coefficients for bulk convection of NaCl in a Brinkman 
medium with charged spherical inclusions: a = 100 nm; I f» 0.951 nm; T = 25° C, 
Di w 1.33 x 10" 9 m^s _1 (Na+); D 2 ~ 2.03 x 10" 9 m 2 s _1 (CP); u* « 5.15 x 10~ 3 m s _1 . 



(e/(kT) 


j—. TJ * / / i m 2 \ 

D eu /(kTa ) 


CjU/(2Ia) 


C jar 






(J = 1,2) 




na — 1 


I = 9.25 x 10~ 6 mol P 1 






— 1 


0. (U X ID 


O ^1 w 1 A — 2 
0. / 1 X 10 


C 1 A w 1 A — 1 

0.14 X 10 


— 2 


i oo % , i aO 

1.32 X 10 


O AO w 1 A — 1 

2.03 x 10 


r* 1 A < / 1A — 1 

5.14 X 10 


—4 


O /I A w 1 aO 

2.44 X 10 


A AC\ w 1A — 1 

4.49 X 10 


r* 1 A . . 1A — 1 

5.14 X 10 


— () 


O 1 A w 1 AO 

3.10 X 10 


O.06 X 10 


C 1/1 w in — 1 

0.14 X 10 


—8 


3.33 x 10° 


6.44 X lO^ 1 


5.14 x 10 _1 


fta = 10 


I = 9.25 x 10" 4 mol P 1 






-1 


4.38 x 10~ 2 


7.76 x 10" 4 


5.14 x 10" 1 


-2 


9.10 x 10~ 2 


2.16 x 10~ 3 


5.14 x 10 _1 


-4 


1.81 x 10 _1 


5.58 x 10~ 3 


5.14 x 10 _1 


-6 


2.19 x lO^ 1 


7.08 x 10~ 3 


5.14 x 10" 1 


-8 


1.96 x lO^ 1 


5.74 x 10~ 3 


5.14 x 10" 1 


/ta = 100 


I = 9.25 x 10~ 2 mol P 1 






-1 


2.55 x 10~ 3 


4.49 x 10~ 6 


5.14 x 10 _1 


-2 


5.33 x 10~ 3 


1.27 x 10~ 5 


5.14 x 10 _1 


-4 


1.14 x 10~ 2 


3.69 x 10~ 5 


5.14 x 10 _1 


-6 


1.64 x 10~ 2 


5.93 x 10^ 5 


5.14 x 10" 1 


-8 


1.70 x 10" 2 


6.14 x 10" 5 


5.14 x 10 _1 



of the electrical power consumption 

V w K°°E 2 AL + 0(0). (3.8) 

Note that the force exerted on the composite 

F = -ApA ps nc{L e /d 4 )A 2 (f>{'S/a 3 )EC E (3.9) 

reflects the pressure required to pump fluid through the network. Therefore, the average 
shear stress r ~ —F/(Ly/A) required to support the composite scales as 

r ~ - V c(L e /d 4 )(A 3/2 /L)(f>{'S/a 3 )EC E . (3.10) 

Since the area A will be set by the flow rate, the length of the composite should be set 
by the maximum allowable shear stress. 

Finally, the pump efficiency, as measured by the ratio of the rate of flow work |QAp| 
to the electrical power consumption V is 

£ ps r ? c(i e /d 4 )[0(3/a 3 )C fJ ] 2 A/(/i oo L). (3.11) 

Because C E depends on na, £, etc., care must be taken in interpreting this equation. 
Nevertheless, geometrical considerations alone clearly favor thin membranes with large 
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cross-sectional area, operating with a low ionic strength (conductivity) and a high inclu- 
sion volume fraction. In practice, an optimal design (with specified flow Q and, perhaps, 
voltage V = Aip) will be constrained by consideration of the mechanical strength of 
the composite (as indicated by Eqns. (|3.9[) and l|3.10[) '). which clearly diminishes with 
decreasing thickness L and increasing area A. 



Now consider the pressure gradient produced by an average electric field with zero av- 
erage flow. This situation may be realized when an electrolyte-saturated composite is 
bounded by a vessel with impenetrable walls. A practical application involves measuring 
the differential (static) pressure Ap to infer the permeability of the polymer gel or, for 
example, the ^-potential of the inclusions. Note that zero average flow does not imply 
stationary fluid at the microscale, because the 'inner' electroosmotic flow around each 
inclusion is balanced by a far-field pressure-driven back-flow, analogous to the situation 
encountered in microelectrophoresis capillaries with blocked ends. 
Setting U = in Eqn. fTGfy gives 



which is termed the incremental pressure gradient. Representative values may be calcu- 
lated by multiplying the incremental pore mobilities in the last column of table ^ and 
plotted in figures |21 EI by their respective values of r\j£ 2 . 

It is interesting to note that when nt, < 1 and, hence, the pore mobility increases 
linearly with £ 2 (see figure |S} , the pressure gradient is independent of the permeabil- 
ity. This is because increasing the permeability increases the electric-ficld-induccd flow, 
and, since the back-flow and accompanying pressure gradient are proportional to each 
other, it follows that the pressure gradient is independent of £ 2 . However, at higher ionic 
strengths, when k£ > 1 (sec figure BJ, the pressure gradient evidently decreases with £ 2 . 
This is because the electric-field-induced flow within the diffuse double layers — where 
resistance to flow is predominantly due to viscous stress — increases more slowly with the 
permeability than the Darcy drag beyond the double layers decreases. 

Note that the incremental pressure gradient reflects the same asymptotic coefficient C E 
as the pore mobility, so it is important to establish whether measuring the electric-ficld- 
induccd pressure gradient offers a significant advantage over measuring the pore mobility. 
Recall, pore mobilities can generate low but measurable velocities in a microchanncl. 
From tableQwith na = 100 and ( = -lkT/e, setting E = 2 V cm -1 and cf> = 10 -2 yields 
P f=a 8.76 kPa cm -1 . Therefore, when L = 5 mm, for example, the pressure differential is 
\Ap\ as 4.38 kPa (static head of 0.43 m of water). Clearly, the pressure gradient induced 
by a relatively weak electric field is sufficient to produce a modest (static) pressure. 



3.3. Incremental pressure gradient 



P = 



c/>{ri/^)(3/a 3 )C E E + 0(<l> 2 ) 



(3.12) 




(3.13) 



where 
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is the sum of incremental microscale contributions (of electromigration, diffusion, and 
convection, respectively) to the average electromigrative flux. 
The ratio of the convective and electromigrative terms is 

A fJ A f,e = — ~ Pe;M, (3-15) 
and the ratio of the convective and diffusive terms is 

Af c /Af d = ~ Pe,, (3.16) 

where the Peclet number Pej = u c K _1 /Z?j is typically very small. The characteristic 
(microscale) velocity u c may be estimated by balancing the O(Ean) electrical force (per 
unit volume) with the 0(u c r]/£ 2 ) Darcy drag force, giving u c ~ Ecrni 2 jr\. With E ~ 
10 2 V cm -1 , a ru 1 ^C cm~ 2 , k" 1 ~ 10 2 nm, I ~ 1 nm, and 77 ~ 10~ 3 kg m _1 s — 1 , 
u c ~ 1 (im s -1 . Further, with a ~ 1 /jm and Dj ~ 10 -9 m 2 s _1 , Pe 3 ~ 10 -4 , indicating 
that diffusion and electromigration (z,j 7^ 0) dominate convection. Clearly, for charged 
species with Zj ~ 1, electromigrative fluxes are comparable to diffusive fluxes. 

Incremental contributions to the ion fluxes are provided in table [3] for the composite 
whose asymptotic coefficients are shown in table ^ These confirm that the convective 
contribution A^ c is small, and that the electromigrative contribution A^ e approaches 
the Maxwell value A^ e = —3/2 for impenetrable (non-conducting) spheres as \Q — > 0. 
Furthermore, the electromigrative (A^ e ) and diffusive (A^ d ) contributions for Na + and 
Cl _ vary significantly with ionic strength. For example, at low electrolyte concentration, 
the diffusive term dominates, enhancing the flux of the counter-ion (Na + ) and attenuating 
that of the co-ion (CI - ). 

3.5. Electrical conductivity 

The electrical conducti vity of colloidal dis persions is well known to reflect the particle 
surface charge density (|Russel et qZ.lll989|) . The conductivity of a composite with im- 
mobilized particles presents a relatively simple problem when the electrolyte ions are 
unhindered by the polymer, because only the electrolyte ions — not the charge on the 
particles (macroions) themselves — contribute to charge transfer. This section establishes 
whether the electrical conductivity is sensitive to the surface charge and, possibly, the 
permeability of the polymer gel. 

From the fluxes in Eqn. 1(3.2(1 . the average current density may be written 

JV N 

I = J2 z i eJ i ~ K °° E + H3/a 3 )E[K°°D E + ^ ZjeDjCf}, (3.17) 
j'=i 3=1 

where 

N n 

K- = ^ef^nf (3.18) 
3=1 

is the conductivity of the electrolyte. The conductivity of the composite is defined as 

K* =I/E = K 00 (1 + 0A K ), (3.19) 

where A K is termed the (dimcnsionless) conductivity increment. 

Equations 1(3. 17(1 - 1(3. 19(1 are equivalent to expressions derived bv lO'"Brienll|l98lh for the 
conductivity of dilute colloidal dispersions with the particles undergoing electrophoresis. 



Table 3. Incremental contributions (see Eqn. 13. 141 ) to the bulk electromigration of NaCl in a Brinkman medium with charged spherical 
inclusions: a = 100 nm; i w 0.951 nm; T = 25°C; £>i w 1.33 x 10~ 9 m 2 s _1 (Na+); D 2 ~ 2.03 x lO^mV 1 (CP). 
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Here, the asymptotic coefficients are different because the particles are stationaryf . When 
the ^-potential is low and, hence, ion fluxes are unperturbed by the surface charge, the 
dipolc strength for non-conducting spheres equals the Maxwell value D E = — (l/2)a 3 , so 
A A — > —3/2 as \C\ — > 0. In general, however, the conductivity increment also reflects the 
charge of the inclusions and the ionic strength. 

Similarly to dispersions, the average convective term (involving C E in Eqn. (|3.2fl 'l 
does not influence the conductivity increment (because of electrical neutrality), and the 
diffusive term (involving C E ) vanishes only when the species have equal mobilities. In 
general, the (microscale) electromigrative and diffusive terms (involving D E and C E ) 
contribute to the average current density. However, because these are influenced by fluid 
motion, the conductivity of a composite is not the same as when particles are fixed in an 
electrolyte without polymer. As expected, because only the electrolyte ions contribute 
to charge transfer, the conductivity of a composite is lower than when particles undergo 
electrophoresis . 

Representative conductivity increments for mobile (solid lines) and stationary (dashed 
lines) particles in a NaCl electrolyte without polymer are compared in figure HJ. At 
low ionic strength (small na), immobilizing the particles decreases the conductivity in- 
crement, because particle migration (electrophoresis) contributes significantly to charge 
transfer. At high ionic strength (na > 10), however, the conductivity increments for 
mobile and fixed particles are (practically) the same. This is because the density of 
charge added by the particles (macroions) is vanishingly small compared to the density 
of bulk electrolyte ions, so the conductivity reflects only the contribution of (dielectric 
and double-layer) polarization to the average electr ic field. 

In contrast to KC1 and HCIO4 electrolytes fsee IO'Brienlll98l]) . the conductivity in- 
crement with NaCl decreases with increasing (^-potential when £ and na are small. This 
may be attributed to the counter- ion (Na + ) having a significantly lower mobility than 
the co- ion (Cl~). For KC1, the mobilities of K + and Cl~ are very similar, yielding a 
monotonically increasing conductivity increment. For HCIO4, however, the counter-ion 
(H + ) has a significantly higher mobility than the co-ion (C10 7), yielding a conductivity 
increment that increases more rapidly (linearly) with |£| fsee lO'Brienlll98l|) . 

The conductivity increment for particles with radius a = 100 nm embedded in a poly- 
mer gel with Brinkman screening length £ sa 3.11 nm (solid lines) is shown in figure 
Similar calculations (not shown) reveal that an order-of-magnitudc increase in the perme- 
ability £ 2 produces almost the same values at all ^-potentials and ionic strengths (values 
of na). This confirms that the average current density is dominated by clcctromigra- 
tion and diffusion. The conductivity increments (solid lines) are compared to values for 
stationary particles in an electrolyte without polymer (dashed lines). These results are 
almost the same at high ionic strength when the diffuse double layer is thin compared 
to the Brinkman screening length, i.e., when k£ 3> 1. Under these conditions, the elec- 
trical force (inside the diffuse double layer) is balanced by viscous stress and, hence, the 
convective flows are similar. When k£ = na(£/a) = 1, for example, na w 100/3.11 « 32 
and, as expected, the limiting behavior (nl ^> 1) at high ionic strength occurs when na 
exceeds this value. 

f When particles undergo electrophoresis in a Newtonian electrolyte, the far-field (solenoidal) 
velocity disturbance is irrotational, decaying as r~ 3 as r — > do. When fixed in electrolyte without 
polymer, however, the far-field disturbance reflects a net force, decaying as r^ 1 as r — * 00. 

X The computations for particles in a pure electrolyte were performed w ith software (called 
MPEK, available from the autho r) based on th e work of iHill et all J2003D . These accurately 
reproduced earlier calculations bv lO'Brienl (|1981ft for KC1 and HCIO4 electrolytes. 
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Figure 6. The (scaled) conductivity increment A K (see Eqn. 13, 191 ) as a function of the (scaled) 
(-potential (e/(kT) for various (scaled) reciprocal double-layer thicknesses na = 1, 2, 10 and 
100 (aqueous NaCl at T — 25° C with a — 100 nm) for particles undergoing electrophoresis (solid 
lines) and stationary particles (dashed lines), both in electrolyte without polymer. 
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Figure 7. The (scaled) conductivity increment A K (see Eqn. 13.191 ') as a function of the 
(scaled) ("-potential Qe/(kT) for various (scaled) reciprocal double-layer thicknesses na = 1, 2, 
10, 100 and 1000: aqueous NaCl at T = 25°C; a = 100 nm; I « 3.11 (solid lines); stationary 
particles in electrolyte without polymer (dashed lines). 



4. Summary 

A rigorous theoretical methodology was presented to calculate steady electrokinctic 
transport of electrolytes in a continuous polymer gel embedded with charged spher- 
ical inclusions. Composites with this microstructure are candidates for enhanced gel- 
electrophoresis, chemical sensing, membrane separation, and, perhaps, electroosmotic 
pumping technologies. This work was also motivated by a desire to interpret experi- 
ments that probe the surface charge of immobilized colloids and the micro-rheology of 
delicate polymer gels. 
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From a numerically exact treatment of electromigration, diffusion, and and convective 
transport past a single inclusion in an unbounded polymer gel, averaged equations de- 
scribing bulk transport properties were derived. The theory was applied to calculate the 
response to a steady electric field with a uniform bulk electrolyte concentration. Note that 
the response to a bulk electrolyte concentration gradient will be treated in a forthcoming 
publication. In this work, electromigration and diffusion were found to be independent of 
convection and, hence, of the polymer-gel (Darcy) permeability. However, the strength 
of clcctroosmotic flow reflects the gel permeability and, to a lesser (but still significant) 
extent, polarization and relaxation by electromigration and diffusion. 

When particles are immobilized in a (neutral) polymer gel, the electrical conductivity is 
independent of the polymer, whereas the pore mobility — similarly to the electrophoretic 
mobility of dispersions — reflects the size and charge of the inclusions and the Darcy 
permeability. Furthermore, when the Debye length is smaller (greater) than the Brinkman 
screening length I (square root of the Darcy permeability), the pore mobility increases 
linearly (quadratically) with I. 

The variation of pore mobility with ^-potential and ionic strength is more compli- 
cated because of the significant influence of polarization and relaxation. Nevertheless, 
the mobility is (approximately) inversely proportional to the inclusion radius, indicating, 
as expected, that the average flow is proportional to the average counter-charge density. 

The Darcy drag of the intervening polymer gel leads to slow flows that will often 
be independent of the differential pressure required to pump fluid through a (modest) 
microfluidic network. Optimal pumping efficiency is favored by thin membranes with 
large cross-section, high inclusion volume fractions, and low electrolyte conductivities. 

The present model assumes that the polymer gel does not influence ion mobilities 
(diffusion coefficients) , and that the volume fraction of the inclusions is small. By analogy 
with Maxwell's well-known theory for conduction in dilute random arrays of spheres, 
the theory advanced in this work may be accurate at moderate volume fractions. Note 
also that the calculations require the bulk electrostatic potential and electrolyte ion 
concentrations to vary slowly in space (and time). Future development of the model will 
accommodate harmonic temporal fluctuations in the applied electric field, permitting the 
interpretation of dielectric relaxation spectroscopy experiments. 
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Appendix A. Derivation of the averaged momentum equation 

This appendix supplements t|2.4l whcre the single-particle microscale problem is adopted 
to quantify the influence of the inclusions on bulk momentum transport. Note that the 
closures in this work neglect hydrodynamic interactions between inclusio ns, wh i ch are 
screened by the intervening Brinkman medium. The reader is referred to iHinchl (119771) 
for details necessary to account for hydrodynamic interactions. Note also that the fields 
obtained from the single-particle problem in this work approximate the conditionally av- 
eraged fields in Hinch's ensemble averaging methodology. Because the microstructure is 
homogeneous, the volume averages below are equivalent to ensemble averages. Integrals 
over the surface or volume of a (spherical) inclusion centered at r = are identified 
by a range of integration that involves the radial coordinate r; otherwise integrals refer 
to "representative control volumes" . Unless stated otherwise, the notation and symbols 
below are the same as in the main text. 

Outside the inclusions, the (inertialess) microscale momentum conservation equation 

is 

= V-T-(r)/l 2 )u-pViP + p f g, (Al) 

where T = —pS + 2r/e is the Newtonian stress tensor, pf is the fluid density, and pj-g is 
a uniform body force (e.g., gravity). 

Similarly, let the microscale momentum conservation equation inside the (rigid) inclu- 
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sions be 

= V-T s +p p g + f 9 , (A2) 

where T s is the stress tensor, and p s is the density. Note that f 9 is a generalized function 
to represent the electrical and mechanical-contact forces acting on the surfaces of the 
inclusions: f 9 — when r ^ a, and 

V- 1 [ f 9 dV = n(f 9 ) (A3) 
Jv 

when nV 1 and V 1 ^ 3 is small compared to the characteristic macroscopic length scale. 

Averaging the momentum conservation equation yields Brinkman's equation for the 
continuous phase (fluid-saturated polymer gel), with additional terms arising from the 
inclusions: 

= V • (T) - ( V /£ 2 )(u) - (pVV) + Pf g 
+ V></ (T s +pS)dV)} + n(f 9 )+cb(p p -p f )g. (A 4) 

Note that (p'Vtp) includes only the counter charge; the influence of the fixed surface 
charge is captured by (f 9 ) , which is the sum of the average electrical (f e ) and mechanical- 
contact (f m ) forces acting on the surfaces of the inclusions (or more generally inside). 
Accordingly, the (averaged) equation of motion for the (immobilized) inclusions is 

n{f d )+n{f g ) + 4>p p g = 0, (A 5) 

where 

(f d ) = (f T-e r dA) (A 6) 

J r—a 

is the average (drag) force exerted by the fluid on the inclusions. 

A useful identity for transforming volume to surface integrals is obtained by differen- 
tiating the product (summation on repeated indices) 

d(akim...Xi)/dxj = Xidakim.../dxj + aki m ..Aj, (A 7) 

where aum... are the components of an arbitrary-order tensor {e.g., stress). Integrating 
this over a volume J dV enclosed by a surface J dA and applying Gauss's integral theorem 
gives 

Ukim...dV = J x l a klm ...h j dA- J x l da k i m .../dx j dV. (A 8) 

Therefore, with T s = T at r — a, and V • T s = —p p g and e = (rigid particles) when 
r < a, it follows that 

(f (T s + P 8)dV) = ( f r2r 1 e-e r dA). (A 9) 

Jr<a J r—a 

In general, body forces and the isotropic stress contribute (J r(f 9 + p p g — 'Vp)dV) to 
the right-hand side of Eqn. IjA 9() . but these integrals vanish if, on average, the (internal) 
body-force distributions and (internal) pressure gradient are even functions of position. 
With a macroscale velocity gradient (not considered in this work), the extensional and 
rotational contributions to e are even functions of position, so (J* _ re ■ e r dA) will be 
linear in (Vu) and, therefore, will contribute to the bulk (deviatoric) stress by modifying 
the effective viscosity rj' (e.g., the well-known Einstein viscosity for dilute, force- free 
suspensions) . 
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Finally, collecting the results from Eqns. IjA 4j) - ljA 9|l gives 

= V • (T) - (ri/e 2 )(u) - ( P Vi>) + (1 - 4>)p f g 
+ V-{n2r 1 {( re ■ e r dA)} - n(f d ) . (A 10) 

J r—a 

The correctness of Eqn. (|A 10|) can be verified, in part, by considering a stationary fixed 
bed of inclusions in the absence of electrical forces. Accordingly, Eqn. IjA 10|) simplifies 
to 

O = -V(p) + (l-0) Pf g-n(f d ), 
and the static cquillibrium of the inclusions requires 

= n(f d )+n(f 9 )+cbp p g. (All) 

Therefore, 

= - V(p) + n(f 9 ) + p f g + <j>{ Pp - Pf )g, (A 12) 

and, hence, in a stationary fluid where V(p) = Pfg, the average force required to immo- 
bilize the inclusions is simply 

(f 9 ) = Wn)( Pf - Pp )g = (4/3)7ra 3 ( P/ - p p )g. (A 13) 



